High Temperature Conductance of the Single Electron Transistor 
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The linear conductance of the single electron transistor is determined in the high temperature 
limit. Electron tunneling is treated nonperturbatively by means of a path integral formulation and 
the conductance is obtained from Kubo's formula. The theoretical predictions are valid for arbitrary 
conductance and found to explain recent experimental data. 
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Considering applications of single electron tunneling, 
one is faced with a dilemma. On the one hand, the con- 
ventional approach Q consider systems with small tun- 
neling conductance Gt *C Gk = e 2 /h implying nearly 
perfect charge quantization. On the other hand, when 
using these devices, e.g. as highly sensitive electrometers 
Q or for thermometry ||, a large current signal is de- 
sirable meaning large tunneling conductance. Since for 
Gt S> Gk charging effects disappear, a compromise must 
be found in practice. This problem has spurred consid- 
erable interest in the precise behavior of single charge 
tunneling devices for large conductance ||,[|. 
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FIG. 1. Circuit diagram of the single electron transistor. 

In this letter, we consider the single electron tran- 
sistor (SET) at high temperatures for arbitrary tunnel- 
ing conductance of the junctions. The SET consists of 
two tunnel junctions in series, with tunneling conduc- 
tance G\ and G 2 and capacitance C\ and C%, respec- 
tively, biased by a voltage source V which may be split 
among the left and right branches into p\V and p 2 V with 
pi + P2 = 1; cf- Fig- §■ The island in between the junc- 
tions is connected via a gate capacitance C g to a con- 
trol voltage U g shifting the electrostatic energy of the 
system continuously. The important energy scale is the 
charging energy E c — e 2 /2C with the island capacitance 
C = G\ + C2 + C g . For weak electron tunneling, E c is 
the energy needed to charge the island with one excess 
electron at vanishing gate voltage U g = 0. Due to the 
periodicity of the Hamiltonian in U g , the conductance is 
a periodic function with period 1 of the dimensionless 



voltage rig = U g C g /e ||. Specifically, we are interested 
in the linear dc conductance for small transport voltage, 
that may be calculated from the Kubo formula 



G = lim - — Im lim 

ui — >0 TlLU iui — ►u?+i<5 



h/3 



dre <1/,T (J (1) (r)/ (2) (0)). (1) 



Since the dc current through both junctions coincides, 
the first current operator JW may be an arbitrary lin- 
ear combination /W = — e 2 I 2 (with £1+62 = 1) 
of the current operators 1\ and I 2 through junctions 1 
and 2, respectively. The relative "— " sign comes from 
the opposite directions of I\ and I 2l which are both 
positive for flux onto the island. The second current 
operator 1^ = k\Ii — K2I2, with the relative weights 
«i = (C2 + piC g )/C and k 2 = (Ci + p 2 C g )/C, is de- 
termined by linear response theory from the coupling of 
the transport voltage V. To evaluate the current-current 
correlator we employ a generating functional 



z[Si,b] = 



{1 r h0 r 
-t / dr iJ-^J^(r) 



, (2) 



where H is the Hamiltonian of the system for V = 0, 
and T T is the time ordering operator for imaginary times 
r. The correlator is then given by a second order varia- 
tional derivative relative to the auxiliary fields £1 and £2 • 
In the phase representation B , we get for the generating 
functional 



with the effective action 

S[p, 6 , 6] = s c [p] + Si [ip, &] + s 2 [ip, U 

The first term on the rhs is the Coulomb action 
V r^ 2(T) 



S c [ip] = dr 



4B„ 



(3) 



(4) 



(5) 



describing the Coulomb charging of the island in presence 
of an applied gate voltage, and the effective tunneling ac- 
tions (i = 1,2) 
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Si[<P,ii]= 9i I drdr' 'a(r - r') 
Jo 

(1 - ie6(r))(l + ieZiWeW-rWV (6) 

describe quasi-particlc tunneling through junctions 1 and 
2, respectively. Here gi — Gi/Gpc is the dimensionless 
conductance of junction i, and the kernel a(r) may be 
written as Fourier series 



1 



(7) 



E = lim — — Im lim . 

UJ — >0 AlU^ ZI/; — 'U + iS Jq 

and 



ai(r)/i(0)) J 



Si 



F=HmJ-Im lim f fr ^M» (14) 
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Since the conductance does not depend on the specific 
choice of the parameters ti and pi , we then find that 



G = E/(G- l + G- 1 ), 



(15) 



where the v n = 2nn/h(3 are Matsubara frequencies. For 
vanishing auxiliary fields the action reduces to the action 
of the single electron box j|] and thus the generating func- 
tional reduces to the box partition function Z = Z[0, 0]. 
Performing the variational derivatives explicitly, we get 
for the correlator 

(WW)) = (/ j (t)/ j -(t')) E ^ + {I,{r)I f {r')) F . (8) 

Since the auxiliary fields are in the argument of an ex- 
ponential, there are two contributions. The first term 
comes from the second order variational derivative of the 
action and reads 



Mt)W)) E = ^9Mt-t')^ J D[ 



e X p|--5M|cos[^(r)-^(rO] I (9) 

with the box action S[tp] = S[<p, £i = 0, £2 = 0]. The sec- 
ond term in Eq. (^) involves a multiplication of two cur- 
rent functionals arising as first order variational deriva- 
tives of the action 
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(WW))* =9i9i'J / D[ 



^ V ^--Sy\\l[ip,T\I[ip,T'l (10) 



with the current functional 
2e f m 

I[<P,t] = — dr'afr - r') sin[c9(r) - u>(/)]. (11) 
K Jo 

Taking into account that {Ij{T)Ij{T r )) E / gj and 
(Ij(r)Iji (t')) f I gjgj' depend only on the dimensionless 
parallel conductance g = gi + gi and thus are indepen- 
dent of the indices j and f, the conductance may be 
written as 

G = eiKxigxE + g\F) - 2e 1 K 2 g\g2F 

+e 2 K 2 (g2E + g 2 F), (12) 

where 



This is a formally exact expression for the linear dc con- 
ductance. 

To proceed, we make explicit the sum over winding 
numbers k of the phase and write the correlator (^|) in 
the form 

tp(K6)=2irk 

O 2 1 00 f 

(7 1 (r)/ 1 (0)) B /. 9l = ^-a(r)- ]T J D[cp] 

k— — oo / n \ n 

expl.-~S[<p}\cos[ip(T)-ip(0)}. (16) 

For given winding number k, the path integral may be 
evaluated approximately by expanding about the classi- 
cal path <p(t) = ip(0) + ffcT. An arbitrary path of wind- 
ing number k is of the form tp(r) — (p(r) + C( T ) with 
C(0) = C(h/3) — 0, and the action can be written in terms 
of the Fourier coefficients £„ = (' n + as 

S[<p + C] = -2irikn g + + S 2 S^ + S*S^ + ... . (17) 

Here 

- 2 k 2 



c(fc) _ 



IE. ' |fc| f' 



(18) 



is the action of the classical path. The second term reads 

00 

6 2 SW = J2\W(C 2 + C), (19) 



n=l 



with the eigenvalues 



X^ = 2 -^+gQ{n-\k\){n-\k\). 
pE c 



(20) 



The term S^S^ is of fourth order in the Fourier com- 
ponents C,n and not given explicitly here. Since the An 
are large for small 0E C , the expansion (|l7| ) about the 
classical path converges rapidly for high temperatures. 

Rewriting the cosine function in Eq. (|l6| ) as a sum of 
exponentials, we get for the correlator the expansion 
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2e 2 1 

(I 1 (r)I 1 (0)) E /g 1 = ^-a(T)^ £ C k e 2 ^ cos(^r) exp 



k=— oo 



1 — cos(^„t) 

An 



1 + 5, 



(fc) 



(21) 



where the coefficients C k read 

= r(i + fc + )r(i + fc_) _ g( « 
fe r 2 (i + fc)r(i + M ) 



(22) 



with fc± = + | ± ±\/4ufc + u 2 and u = g(5E c /2'n: 2 . The 
dominant corrections to the semiclassical approximation 
are described by 



,(fc) .(fc) 
A m A, 



6t(i/jt) - 2a(i/„+fe) 



-2a(i/„_ fc ) + a(v m+n+k ) + &(v m+n -k) ■ (23) 



(24) 



with the same correction (^3|). The expansions ( pl| ) and 
( p4| ) proceed in powers of /3£cj however, terms involv- 
ing u = gl3E c /2ir 2 are kept to a/Z orders. This ensures 
a meaningful result in the limit of moderately high tem- 
peratures also for large conductance g. 

From Eq. (Ell) one obtains for the Fourier coefficients 



dre i ^{h{T)I 1 {Q)) E /g 1 



(25) 



The corresponding expansion of the partition function Z 
reads 



the high temperature expansion 



exp 



n— 1 



+ 2 — A w A w ha(,^+fep 4 +. 

771, n^O m 71 



rn/0 



(26) 



When these coefficients are analytically continued in 
the complex v plane, they are analytic on each half plane 
Rez;>0 with a cut along the imaginary axis [||. The 
representation of E(v) as a sum over winding numbers 
k shifts this cut to Re^ = k for the fcth term of the 
sum. Thus, in the phase representation, only the full 
sum shows the analytic properties underlying the con- 
ductance formula (Q) . To deal with this problem we for- 
mally change to the charge representation, perform the 
analytic continuation and the ui —> limit there, and 
then go back to the phase representation. This way the 
high temperature expansion of the conductance (|l^) may 
be evaluated to read 

G - [1/(G? + G^Z- 1 exp {-2[ 7 + + u)\/g} 

{1-^'(1+U)(/?£ C /7T 2 ) 

+ \ga{u) + t(u)} {(3E c /2tt 2 ) 2 + 0{(3E c f). (27) 



3(u,v) 



j/jl+u) 
u hx{v) 

2(1 - v u 



1 



1 - 2v v 
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The dependence on u = gl3E c /2ir 2 is given in terms of 
the digamma function ip(z) and two auxiliary functions 



a{u) 



and 



t(u) = —ip(l + u) - 
bu 



7 + V(l + u) - uip'(l + u) 

\ 2v(l-v u )<j)(v,l,l + u) 
dv 

(1 — VjU 

[</>(! +n)+ 7 ] 2 



(28) 



+ / dvS(u, v) 



3 7 ^'(1 + u) + i>(l + u)[\ + 2^/(1 + u)] 



with Lerch's transcendent 



u 

(z,u,v) 



(29) 



and 



(1 — v)u 
[ln(l -v) + v(j)(v, 1, 1 



ln(«)^(-» 1.1 + «) + #-, 2,1 
v v 



■2v<j>(v,2,l + u) 



•«)] 



21n(w)ln(l - v) + -hi 2 (v) + 3Li 2 (l - v) 



(30) 



The high temperature expansion of Z is straightfor- 
ward and reads 

Z = l + ga{u)(PE c /27T 2 ) 2 + 0{(3E C ) 3 

+2Ci cos(27rn g ) [l + 0(f3E c ) 2 ] , (31) 

which combines with Eq. J27|) to yield an analytical ex- 



pression for the high temperature conduction of a SET 
valid for arbitrary tunneling conductance. 

In the region of weak tunneling, g < 1, the quantity u 
becomes small at high temperatures and we may replace 
<j{u) and t(u) by 



3 



o-(O) = 3C(3), r(0) = 7r7l0. 



(32) 



This gives for the conductance of a weakly conducting 
SET 



G = 



1 



cr 1 



1 - 



1 

15 



3 

3C(3) 



2tt 4 



(33) 



in accordance with earlier work p0[ . In the region of 
strong tunneling the quantity u is typically large even 
for the highest temperatures explored experimentally and 
the full expression (pTj), (|3l|) must be used. 
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FIG. 2. Maximum and minimum linear conductance in de- 
pendence on dimensionless temperature for two dimensionless 
parallel conductances a = 2.5 and 7.3 compared with experi- 
mental data (symbols) by Joyez et al. JlT| , 

We have compared our findings with recent experimen- 
tal data by Joyez et al. jll|] for transistors with g = 0.6, 
2.5 and 7.3. As seen from Fig. the theory describes the 
high temperature behavior of all junctions (results for 
g = 0.6 are not shown) down to temperatures where the 
current starts to modulate with the gate voltage. The 
parameters have not been adjusted to improve the fit 



but coincide with the values given in @. The small 
deviations between theory and experiment for g = 7.3 
near (3E C = 1 may arise from experimental uncertain- 
ties in (3E C . We mention that the temperature depen- 
dence of the conductance of the highly conducting SET 
(g = 7.3) is not within reach of previous theoretical pre- 
dictions. The results obtained thus present substantial 
progress and should be useful for experimental studies of 
even larger tunneling conductances since the predictions 
remain valid for arbitrary values of g. 
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